Cysteine desulfurase (IscS)-mediated fine-tuning of bioenergetics and SUF expression prevents Mycobacterium tuberculosis hypervirulence.

Iron-sulfur (Fe-S) biogenesis requires multiprotein assembly systems, SUF and ISC, in most prokaryotes. M. tuberculosis (Mtb) encodes a complete SUF system, the depletion of which was bactericidal. The ISC operon is truncated to a single gene iscS (cysteine desulfurase), whose function remains uncertain. Here, we show that MtbΔiscS is bioenergetically deficient and hypersensitive to oxidative stress, antibiotics, and hypoxia. MtbΔiscS resisted killing by nitric oxide (NO). RNA sequencing indicates that IscS is important for expressing regulons of DosR and Fe-S-containing transcription factors, WhiB3 and SufR. Unlike wild-type Mtb, MtbΔiscS could not enter a stable persistent state, continued replicating in mice, and showed hypervirulence. The suf operon was overexpressed in MtbΔiscS during infection in a NO-dependent manner. Suppressing suf expression in MtbΔiscS either by CRISPR interference or upon infection in inducible NO-deficient mice arrests hypervirulence. Together, Mtb redesigned the ISC system to "fine-tune" the expression of SUF machinery for establishing persistence without causing detrimental disease in the host.

Understanding the Genome-Wide Transcription Response To Various cAMP Levels in Bacteria Using Phenomenological Models.

Attempts to understand gene regulation by global transcription factors have largely been limited to expression studies under binary conditions of presence and absence of the transcription factor. Studies addressing genome-wide transcriptional responses to changing transcription factor concentration at high resolution are lacking. Here, we create a data set containing the entire Escherichia coli transcriptome in Luria-Bertani (LB) broth as it responds to 10 different cAMP concentrations spanning the biological range. We use the Hill's model to accurately summarize individual gene responses into three intuitively understandable parameters, Emax, n, and k, reflecting the sensitivity, nonlinearity, and midpoint of the dynamic range. Our data show that most cAMP-regulated genes have an n of >2, with their k values centered around the wild-type concentration of cAMP. Additionally, cAMP receptor protein (CRP) affinity to a promoter is correlated with Emax but not k, hinting that a high-affinity CRP promoter need not ensure transcriptional activation at lower cAMP concentrations and instead affects the magnitude of the response. Finally, genes belonging to different functional classes are tuned to have different k, n, and Emax values. We demonstrate that phenomenological models are a better alternative for studying gene expression trends than classical clustering methods, with the phenomenological constants providing greater insights into how genes are tuned in a regulatory network. IMPORTANCE Different genes may follow different trends in response to various transcription factor concentrations. In this study, we ask two questions: (i) what are the trends that different genes follow in response to changing transcription factor concentrations and (ii) what methods can be used to extract information from the gene trends so obtained. We demonstrate a method to analyze transcription factor concentration-dependent genome-wide expression data using phenomenological models. Conventional clustering methods and principal-component analysis (PCA) can be used to summarize trends in data but have limited interpretability. The use of phenomenological models greatly enhances the interpretability and thus utility of conventional clustering. Transformation of dose-response data into phenomenological constants opens up avenues to ask and answer many different kinds of question. We show that the phenomenological constants obtained from the model fits can be used to generate insights about network topology and allows integration of other experimental data such as chromatin immunoprecipitation sequencing (ChIP-seq) to understand the system in greater detail.

Emergence of networks of shared restriction-modification systems in phage-bacteria ecosystems.

Restriction-modification (RM) systems are the most ubiquitous bacterial defence systems against bacteriophages. Using genome sequence data, we showed that RM systems are often shared among bacterial strains in a structured way. Examining the network of interconnections between bacterial strains within genera, we found that many strains share more RM systems than expected compared with a suitable null model. We also found that many genera have a larger than expected number of bacterial strains with unique RM systems. We used population dynamics models of closed and open phage-bacteria ecosystems to qualitatively understand the selection pressures that could lead to such network structures with enhanced overlap or uniqueness. In our models, we found that the phages impose a selection pressure that favours bacteria with greater number of RM systems, and higher overlap of RM systems with other strains, but in bacteria-dominated states, this is opposed by the increased cost-to-growth rate of these bacteria. Similar to what we observed in the genome data, we found that two distinct bacterial strategies emerge - strains either have a greater overlap than expected, or, at the other extreme, have unique RM systems. The former strategy appears to dominate when the repertoire of available RM systems is smaller but the average number of RM systems per strain is larger.

Replication-Dependent Organization Constrains Positioning of Long DNA Repeats in Bacterial Genomes.

Bacterial genome organization is primarily driven by chromosomal replication from a single origin of replication. However, chromosomal rearrangements, which can disrupt such organization, are inevitable in nature. Long DNA repeats are major players mediating rearrangements, large and small, via homologous recombination. Since changes to genome organization affect bacterial fitness-and more so in fast-growing than slow-growing bacteria-and are under selection, it is reasonable to expect that genomic positioning of long DNA repeats is also under selection. To test this, we identified identical DNA repeats of at least 100 base pairs across ∼6,000 bacterial genomes and compared their distribution in fast- and slow-growing bacteria. We found that long identical DNA repeats are distributed in a non-random manner across bacterial genomes. Their distribution differs in the overall number, orientation, and proximity to the origin of replication, between fast- and slow-growing bacteria. We show that their positioning-which might arise from a combination of the processes that produce repeats and selection on rearrangements that recombination between repeat elements might cause-permits less disruption to the replication-dependent genome organization of bacteria compared with random suggesting it as a major constraint to positioning of long DNA repeats.

Evolutionary and Comparative Analysis of Bacterial Nonhomologous End Joining Repair.

DNA double-strand breaks (DSBs) are a threat to genome stability. In all domains of life, DSBs are faithfully fixed via homologous recombination. Recombination requires the presence of an uncut copy of duplex DNA which is used as a template for repair. Alternatively, in the absence of a template, cells utilize error-prone nonhomologous end joining (NHEJ). Although ubiquitously found in eukaryotes, NHEJ is not universally present in bacteria. It is unclear as to why many prokaryotes lack this pathway. Toward understanding what could have led to the current distribution of bacterial NHEJ, we carried out comparative genomics and phylogenetic analysis across ∼6,000 genomes. Our results show that this pathway is sporadically distributed across the phylogeny. Ancestral reconstruction further suggests that NHEJ was absent in the eubacterial ancestor and can be acquired via specific routes. Integrating NHEJ occurrence data for archaea, we also find evidence for extensive horizontal exchange of NHEJ genes between the two kingdoms as well as across bacterial clades. The pattern of occurrence in bacteria is consistent with correlated evolution of NHEJ with key genome characteristics of genome size and growth rate; NHEJ presence is associated with large genome sizes and/or slow growth rates, with the former being the dominant correlate. Given the central role these traits play in determining the ability to carry out recombination, it is possible that the evolutionary history of bacterial NHEJ may have been shaped by requirement for efficient DSB repair.

Whole-Genome Sequencing of Sphingobium sp. Strain RSMS, a Highly Efficient Tributyl Phosphate-Degrading Bacterium.

Sphingobium sp. strain RSMS was described earlier as an efficient degrader of tributyl phosphate, an organic pollutant. This report describes the generation and annotation of the genome sequence of Sphingobium sp. strain RSMS, which will facilitate future studies to identify genetic elements responsible for the degradation of tributyl phosphate.

Dynamics of genetic variation in transcription factors and its implications for the evolution of regulatory networks in Bacteria.

The evolution of regulatory networks in Bacteria has largely been explained at macroevolutionary scales through lateral gene transfer and gene duplication. Transcription factors (TF) have been found to be less conserved across species than their target genes (TG). This would be expected if TFs accumulate mutations faster than TGs. This hypothesis is supported by several lab evolution studies which found TFs, especially global regulators, to be frequently mutated. Despite these studies, the contribution of point mutations in TFs to the evolution of regulatory network is poorly understood. We tested if TFs show greater genetic variation than their TGs using whole-genome sequencing data from a large collection of Escherichia coli isolates. TFs were less diverse than their TGs across natural isolates, with TFs of large regulons being more conserved. In contrast, TFs showed higher mutation frequency in adaptive laboratory evolution experiments. However, over long-term laboratory evolution spanning 60 000 generations, mutation frequency in TFs gradually declined after a rapid initial burst. Extrapolating the dynamics of genetic variation from long-term laboratory evolution to natural populations, we propose that point mutations, conferring large-scale gene expression changes, may drive the early stages of adaptation but gene regulation is subjected to stronger purifying selection post adaptation.

Laboratory Evolution Experiments Help Identify a Predominant Region of Constitutive Stable DNA Replication Initiation.

The bacterium Escherichia coli can initiate replication in the absence of the replication initiator protein DnaA and/or the canonical origin of replication oriC in a ΔrnhA background. This phenomenon, which can be primed by R-loops, is called constitutive stable DNA replication (cSDR). Whether DNA replication during cSDR initiates in a stochastic manner through the length of the chromosome or at specific sites and how E. coli can find adaptations to loss of fitness caused by cSDR remain inadequately answered. We use laboratory evolution experiments of ΔrnhA-ΔdnaA strains followed by deep sequencing to show that DNA replication preferentially initiates within a broad region located ∼0.4 to 0.7 Mb clockwise of oriC. This region includes many bisulfite-sensitive sites, which have been previously defined as R-loop-forming regions, and includes a site containing sequence motifs that favor R-loop formation. Initiation from this region would result in head-on replication-transcription conflicts at rRNA loci. Inversions of these rRNA loci, which can partly resolve these conflicts, help the bacterium suppress the fitness defects of cSDR. These inversions partially restore the gene expression changes brought about by cSDR. The inversion, however, increases the possibility of conflicts at essential mRNA genes, which would utilize only a minuscule fraction of RNA polymerase molecules, most of which transcribe rRNA genes. Whether subsequent adaptive strategies would attempt to resolve these conflicts remains an open question.IMPORTANCE The bacterium E. coli can replicate its DNA even in the absence of the molecules that are required for canonical replication initiation. This often requires the formation of RNA-DNA hybrid structures and is referred to as constitutive stable DNA replication (cSDR). Where on the chromosome does cSDR initiate? We answer this question using laboratory evolution experiments and genomics and show that selection favors cSDR initiation predominantly at a region ∼0.6 Mb clockwise of oriC. Initiation from this site will result in more head-on collisions of DNA polymerase with RNA polymerase operating on rRNA loci. The bacterium adapts to this problem by inverting a region of the genome including several rRNA loci such that head-on collisions between the two polymerases are minimized. Understanding such evolutionary strategies in the context of cSDR can provide insights into the potential causes of resistance to antibiotics that target initiation of DNA replication.

Early fate of exogenous promoters in E. coli.

Gene gain by horizontal gene transfer is a major pathway of genome innovation in bacteria. The current view posits that acquired genes initially need to be silenced and that a bacterial chromatin protein, H-NS, plays a role in this silencing. However, we lack direct observation of the early fate of a horizontally transferred gene to prove this theory. We combine sequencing, flow cytometry and sorting, followed by microscopy to monitor gene expression and its variability after large-scale random insertions of a reporter gene in a population of Escherichia coli bacteria. We find that inserted promoters have a wide range of gene-expression variability related to their location. We find that high-expression clones carry insertions that are not correlated with H-NS binding. Conversely, binding of H-NS correlates with silencing. Finally, while most promoters show a common level of extrinsic noise, some insertions show higher noise levels. Analysis of these high-noise clones supports a scenario of switching due to transcriptional interference from divergent ribosomal promoters. Altogether, our findings point to evolutionary pathways where newly-acquired genes are not necessarily silenced, but may immediately explore a wide range of expression levels to probe the optimal ones.

Targeting redox heterogeneity to counteract drug tolerance in replicating Mycobacterium tuberculosis.

The capacity of Mycobacterium tuberculosis (Mtb) to tolerate multiple antibiotics represents a major problem in tuberculosis (TB) management. Heterogeneity in Mtb populations is one of the factors that drives antibiotic tolerance during infection. However, the mechanisms underpinning this variation in bacterial population remain poorly understood. Here, we show that phagosomal acidification alters the redox physiology of Mtb to generate a population of replicating bacteria that display drug tolerance during infection. RNA sequencing of this redox-altered population revealed the involvement of iron-sulfur (Fe-S) cluster biogenesis, hydrogen sulfide (H2S) gas, and drug efflux pumps in antibiotic tolerance. The fraction of the pH- and redox-dependent tolerant population increased when Mtb infected macrophages with actively replicating HIV-1, suggesting that redox heterogeneity could contribute to high rates of TB therapy failure during HIV-TB coinfection. Pharmacological inhibition of phagosomal acidification by the antimalarial drug chloroquine (CQ) eradicated drug-tolerant Mtb, ameliorated lung pathology, and reduced postchemotherapeutic relapse in in vivo models. The pharmacological profile of CQ (C max and AUClast) exhibited no major drug-drug interaction when coadministered with first line anti-TB drugs in mice. Our data establish a link between phagosomal pH, redox metabolism, and drug tolerance in replicating Mtb and suggest repositioning of CQ to shorten TB therapy and achieve a relapse-free cure.

Regulation of Global Transcription in Escherichia coli by Rsd and 6S RNA.

In Escherichia coli, the sigma factor σ70 directs RNA polymerase to transcribe growth-related genes, while σ38 directs transcription of stress response genes during stationary phase. Two molecules hypothesized to regulate RNA polymerase are the protein Rsd, which binds to σ70, and the non-coding 6S RNA which binds to the RNA polymerase-σ70 holoenzyme. Despite multiple studies, the functions of Rsd and 6S RNA remain controversial. Here we use RNA-Seq in five phases of growth to elucidate their function on a genome-wide scale. We show that Rsd and 6S RNA facilitate σ38 activity throughout bacterial growth, while 6S RNA also regulates widely different genes depending upon growth phase. We discover novel interactions between 6S RNA and Rsd and show widespread expression changes in a strain lacking both regulators. Finally, we present a mathematical model of transcription which highlights the crosstalk between Rsd and 6S RNA as a crucial factor in controlling sigma factor competition and global gene expression.

Using stochastic cell division and death to probe minimal units of cellular replication.

The invariant cell initiation mass measured in bacterial growth experiments has been interpreted as a minimal unit of cellular replication. Here we argue that the existence of such minimal units induces a coupling between the rates of stochastic cell division and death. To probe this coupling we tracked live and dead cells in Escherichia coli populations treated with a ribosome-targeting antibiotic. We find that the growth exponent from macroscopic cell growth or decay measurements can be represented as the difference of microscopic first-order cell division and death rates. The boundary between cell growth and decay, at which the number of live cells remains constant over time, occurs at the minimal inhibitory concentration (MIC) of the antibiotic. This state appears macroscopically static but is microscopically dynamic: division and death rates exactly cancel at MIC but each is remarkably high, reaching 60% of the antibiotic-free division rate. A stochastic model of cells as collections of minimal replicating units we term 'widgets' reproduces both steady-state and transient features of our experiments. Sub-cellular fluctuations of widget numbers stochastically drive each new daughter cell to one of two alternate fates, division or death. First-order division or death rates emerge as eigenvalues of a stationary Markov process, and can be expressed in terms of the widget's molecular properties. High division and death rates at MIC arise due to low mean and high relative fluctuations of widget number. Isolating cells at the threshold of irreversible death might allow molecular characterization of this minimal replication unit.

Repression of YdaS Toxin Is Mediated by Transcriptional Repressor RacR in the Cryptic rac Prophage of Escherichia coli K-12.

Horizontal gene transfer is a major driving force behind the genomic diversity seen in prokaryotes. The cryptic rac prophage in Escherichia coli K-12 carries the gene for a putative transcription factor RacR, whose deletion is lethal. We have shown that the essentiality of racR in E. coli K-12 is attributed to its role in transcriptionally repressing toxin gene(s) called ydaS and ydaT, which are adjacent to and coded divergently to racR. IMPORTANCE Transcription factors in the bacterium E. coli are rarely essential, and when they are essential, they are largely toxin-antitoxin systems. While studying transcription factors encoded in horizontally acquired regions in E. coli, we realized that the protein RacR, a putative transcription factor encoded by a gene on the rac prophage, is an essential protein. Here, using genetics, biochemistry, and bioinformatics, we show that its essentiality derives from its role as a transcriptional repressor of the ydaS and ydaT genes, whose products are toxic to the cell. Unlike type II toxin-antitoxin systems in which transcriptional regulation involves complexes of the toxin and antitoxin, repression by RacR is sufficient to keep ydaS transcriptionally silent.

CRISPR-Cas-Mediated Gene Silencing Reveals RacR To Be a Negative Regulator of YdaS and YdaT Toxins in Escherichia coli K-12.

Bacterial genomes are rich in horizontally acquired prophages. racR is an essential gene located in the rac prophage that is resident in many Escherichia coli genomes. Employing a clustered regularly interspaced short palindromic repeat (CRISPR)-Cas-based gene silencing approach, we show that RacR is a negative regulator of the divergently transcribed and adjacent ydaS-ydaT operon in Escherichia coli K-12. Overexpression of YdaS and YdaT due to RacR depletion leads to cell division defects and decrease in survival. We further show that both YdaS and YdaT can act independently as toxins and that RacR serves to counteract the toxicity by tightly downregulating the expression of these toxins. IMPORTANCEracR is an essential gene and one of the many poorly studied genes found on the rac prophage element that is present in many Escherichia coli genomes. Employing a CRISPR-based approach, we have silenced racR expression to various levels and elucidated its physiological consequences. We show that the downregulation of racR leads to upregulation of the adjacent ydaS-ydaT operon. Both YdaS and YdaT act as toxins by perturbing the cell division resulting in enhanced cell killing. This work establishes a physiological role for RacR, which is to keep the toxic effects of YdaS and YdaT in check and promote cell survival. We, thus, provide a rationale for the essentiality of racR in Escherichia coli K-12 strains.

Modulation of Global Transcriptional Regulatory Networks as a Strategy for Increasing Kanamycin Resistance of the Translational Elongation Factor-G Mutants in Escherichia coli.

Evolve and resequence experiments have provided us a tool to understand bacterial adaptation to antibiotics. In our previous work, we used short-term evolution to isolate mutants resistant to the ribosome targeting antibiotic kanamycin, and reported that Escherichia coli develops low cost resistance to kanamycin via different point mutations in the translation Elongation Factor-G (EF-G). Furthermore, we had shown that the resistance of EF-G mutants could be increased by second site mutations in the genes rpoD/cpxA/topA/cyaA Mutations in three of these genes had been discovered in earlier screens for aminoglycoside resistance. In this work, we expand our understanding of these second site mutations, the goal being to understand how these mutations affect the activities of the mutated gene products to confer resistance. We show that the mutation in cpxA most likely results in an active Cpx stress response. Further evolution of an EF-G mutant in a higher concentration of kanamycin than what was used in our previous experiments identified the cpxA locus as a primary target for a significant increase in resistance. The mutation in cyaA results in a loss of catalytic activity and probably results in resistance via altered CRP function. Despite a reduction in cAMP levels, the CyaAN600Y mutant has a transcriptome indicative of increased CRP activity, pointing to an unknown role for CyaA and / or cAMP in gene expression. From the transcriptomes of double and single mutants, we describe the epistasis between the mutation in EF-G and these second site mutations. We show that the large scale transcriptomic changes in the topoisomerase I (FusAA608E-TopAS180L) mutant likely result from increased negative supercoiling in the cell. Finally, genes with known roles in aminoglycoside resistance were present among the misregulated genes in the mutants.

Nucleoid-Associated Proteins: Genome Level Occupancy and Expression Analysis.

The advent of Chromatin Immunoprecipitation sequencing (ChIP-Seq) has allowed the identification of genomic regions bound by a DNA binding protein in-vivo on a genome-wide scale. The impact of the DNA binding protein on gene expression can be addressed using transcriptome experiments in appropriate genetic settings. Overlaying the above two sources of data enables us to dissect the direct and indirect effects of a DNA binding protein on gene expression. Application of these techniques to Nucleoid Associated Proteins (NAPs) and Global Transcription Factors (GTFs) has underscored the complex relationship between DNA-protein interactions and gene expression change, highlighting the role of combinatorial control. Here, we demonstrate the usage of ChIP-Seq to infer binding properties and transcriptional effects of NAPs such as Fis and HNS, and the GTF CRP in the model organism Escherichia coli K12 MG1655 (E. coli).

Genomewide Mutational Diversity in Escherichia coli Population Evolving in Prolonged Stationary Phase.

Prolonged stationary phase is an approximation of natural environments presenting a range of stresses. Survival in prolonged stationary phase requires alternative metabolic pathways for survival. This study describes the repertoire of mutations accumulating in starving Escherichia coli populations in lysogeny broth. A wide range of mutations accumulates over the course of 1 month in stationary phase. Single nucleotide polymorphisms (SNPs) constitute 64% of all mutations. A majority of these mutations are nonsynonymous and are located at conserved loci. There is an increase in genetic diversity in the evolving populations over time. Computer simulations of evolution in stationary phase suggest that the maximum frequency of mutations observed in our experimental populations cannot be explained by neutral drift. Moreover, there is frequent genetic parallelism across populations, suggesting that these mutations are under positive selection. Finally, functional analysis of mutations suggests that regulatory mutations are frequent targets of selection. IMPORTANCE Prolonged stationary phase in bacteria, contrary to its name, is highly dynamic, with extreme nutrient limitation as a predominant stress. Stationary-phase cultures adapt by rapidly selecting a mutation(s) that confers a growth advantage in stationary phase (GASP). The phenotypic diversity of starving E. coli populations has been studied in detail; however, only a few mutations that accumulate in prolonged stationary phase have been described. This study documented the spectrum of mutations appearing in Escherichia coli during 28 days of prolonged starvation. The genetic diversity of the population increases over time in stationary phase to an extent that cannot be explained by random, neutral drift. This suggests that prolonged stationary phase offers a great model system to study adaptive evolution by natural selection.

Comparative Genomics of Interreplichore Translocations in Bacteria: A Measure of Chromosome Topology?

Genomes evolve not only in base sequence but also in terms of their architecture, defined by gene organization and chromosome topology. Whereas genome sequence data inform us about the changes in base sequences for a large variety of organisms, the study of chromosome topology is restricted to a few model organisms studied using microscopy and chromosome conformation capture techniques. Here, we exploit whole genome sequence data to study the link between gene organization and chromosome topology in bacteria. Using comparative genomics across ∼250 pairs of closely related bacteria we show that: (a) many organisms show a high degree of interreplichore translocations throughout the chromosome and not limited to the inversion-prone terminus (ter) or the origin of replication (oriC); (b) translocation maps may reflect chromosome topologies; and (c) symmetric interreplichore translocations do not disrupt the distance of a gene from oriC or affect gene expression states or strand biases in gene densities. In summary, we suggest that translocation maps might be a first line in defining a gross chromosome topology given a pair of closely related genome sequences.

The genome-scale interplay amongst xenogene silencing, stress response and chromosome architecture in Escherichia coli.

The gene expression state of exponentially growing Escherichia coli cells is manifested by high expression of essential and growth-associated genes and low levels of stress-related and horizontally acquired genes. An important player in maintaining this homeostasis is the H-NS-StpA gene silencing system. A Δhns-stpA deletion mutant results in high expression of otherwise-silent horizontally acquired genes, many located in the terminus-half of the chromosome, and an indirect downregulation of many highly expressed genes. The Δhns-stpA double mutant displays slow growth. Using laboratory evolution we address the evolutionary strategies that E. coli would adopt to redress this gene expression imbalance. We show that two global gene regulatory mutations-(i) point mutations inactivating the stress-responsive sigma factor RpoS or σ38 and (ii) an amplification of ∼40% of the chromosome centred around the origin of replication-converge in partially reversing the global gene expression imbalance caused by Δhns-stpA. Transcriptome data of these mutants further show a three-way link amongst the global gene regulatory networks of H-NS and σ38, as well as chromosome architecture. Increasing gene expression around the terminus of replication results in a decrease in the expression of genes around the origin and vice versa; this appears to be a persistent phenomenon observed as an association across ∼300 publicly-available gene expression data sets for E. coli. These global suppressor effects are transient and rapidly give way to more specific mutations, whose roles in reversing the growth defect of H-NS mutations remain to be understood.